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TO THE MEMORY OF MARTIN GARDNER (OCTOBER 21, 1914 - MAY 22, 2010). 

July 5, 2011 

Abstract 

Dismal arithmetic is just like the arithmetic you learned in school, only simpler: there there 
are no carries, when you add digits you just take the largest, and when you multiply digits 
you take the smallest. This paper studies basic number theory in this world, including 
analogues of the primes, number of divisors, sum of divisors, and the partition function. 
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To remedy the dismal state of arithmetic skills possessed by today's children, we propose a 
"dismal arithmetic" that will be easier to learn than the usual version. It is easier because 
there are no carry digits and there is no need to add or multiply digits, or to do anything 
harder than comparing. In dismal arithmetic, for each pair of digits, 

to Add, take the lArger, but 

to Multiply, take the sMaller. 
That's it! For example: 2 + 5 = 5, 2x5 = 2. 

Addition or multiplication of larger numbers uses the same rules, always with the proviso 
that there are no carries. For example, the dismal sum of 169 and 248 is 269 and their dismal 
product is 12468 (Figure 1). 
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Fig. 1(a) Dismal addition. 



Fig. 1(b) Dismal multiplication. 

One might expect that nothing interesting could arise from such simple rules. However, 
developing the dismal analogue of ordinary elementary number theory will lead us to some 
surprisingly difficult questions. 

Here are a few dismal analogues of standard sequences. The "even" numbers, 2xn, are 

0, 1, 2, 2, 2, 2, 2, 2, 2, 2, 10, 11, 12, 12, 12, 12, 12, 12, 12, 12, 20, 21, 22, 22, ■ ■ • (1) 

(entry A171818 in [17]). Note that n + n (which is simply n) is a different sequence. For 
another, less obvious, analogue of the even numbers, see (13) in §3. The squares, nx n, are 

0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 100, 111, 112, 113, 114, 115, 116, 117, 118, 119, 200, • • ■ (2) 

(A087019), 1 the dismal triangular numbers, + 1 + 2 . . . + n, are 

0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 29, 29, 29, 29, 29, . . . (3) 

(A087052), and the dismal factorials, lx2x •■• xn, n>l, are 

1, 1, 1, 1, 1, 1, 1, 1, 1, 10, 110, 1110, 11110, 111110, 1111110, 11111110, 111111110, 
1111111110, 11111111110, 111111111100, 1111111111100, 11111111111100, ... (4) 

(A189788). 

A formal definition of dismal arithmetic is given in §2, valid for any base b, not just 
base 10, and it shown there that the commutative, associative, and distributive laws hold 
(Theorem 1). In that section we also introduce the notion of a "digit map," in order to 
study how changing individual digits in a dismal calculation affects the answer (Theorem 3, 
Corollary 4). 

The dismal primes are the subject of §3. A necessary condition for a number to be a 
prime is that it contain a digit equal to b — 1. The data suggest that if k is large, almost all 
numbers of length k containing b — 1 as a digit and not ending with zero are prime, and so 
the number of primes of length k appears to approach {b — l) 2 b k ~ 2 as k — > oo (Conjecture 1). 
In any case, any number with a digit equal to b — 1 is a product of primes (Theorem 7), and 
every number can be written as r times a product of primes, for some r G {0, 1, . . . , b — 1} 



1 See also the sums of two squares, A171120. The numbers 10, 11, . . . , 99 are not the sum of any number 
of squares, so there is no dismal analogue of the four-squares theorem. 



(Corollary 8). These factorizations are in general not unique. There is a useful process using 
digit maps for "promoting" a prime from a lower base to a higher base, which enables us to 
replace the list of all primes by a shorter list of prime "templates" (Table 3). 

Dismal squares are briefly discussed in §4. 

In §5 we investigate the different ways to order the dismal numbers, and in particular 
the partially ordered set defined by the divisibility relation (see Table 1). We will see that 
greatest common divisors and least common multiples need not exist, so this poset fails to 
be a lattice. On the other hand, we do have the notion of "relatively prime" and we can 
define an analogue of the Euler totient function. 

In §6 we study the number-of-divisors function db(n), and investigate which numbers 
have the most divisors. It appears that in any base b > 3, the number n = (b k — l)/(fe— 1) = 
111 ... 1 It has more divisors than any other number of length k. The binary case is slightly 
different. Here it appears that among all fc-digit numbers n, the maximal value of d 2 {n) 
occurs at n = 2 k — 2 = 111 . . . IOJ2, and this is the unique maximum for n 7^ 2,4. Among 
all odd /c-digit numbers n, ^(w) has a unique maximum at n = 2 k — 1 = 111... Ill I2, 
and if k > 3 and k 7^ 5, the next largest value occurs at n = 2 k — 3 = 111 . . . 101 12, its 
reversal 2 k — 2 k ~ 2 — 1 = 101 . . . 1 11 1 2 , and possibly other values of n (see Conjectures 2-4). 
Although we cannot prove these conjectures, we are able to determine the exact values of 
c4(lll . . . lll|fe) and o?2(lll • • • IOII2) (Theorem 13, which extends earlier work of Richard 
Schroeppel and the second author, and Theorem 14). 

The sequence of the number of divisors of 11 . . . 1 1 1 2 (with k l's) turns out to arise in a 
variety of different problems, involving compositions, trees, polyominoes, Dyck paths, etc.— 
see Remark (iii) following Theorem 12. The initial terms can be seen in Table 8. This 
sequence appears in two entries in [17], A007059 and A079500, and is the subject of a survey 
article by Frosini and Rinaldi [6]. The asymptotic behavior of this sequence was determined 
by Kemp [10] and by Knopfmacher and Robbins [12], the latter using the method of Mellin 
transforms — see (23). This is an example of an asymptotic expansion where the leading 
term has an oscillating component which, though small, does not go to zero. It is amusing 
to note that one of the first problems in which the asymptotic behavior was shown to involve 
a nonvanishing oscillating term was the analysis of the average number of carries when two 
/c-digit numbers are added (Knuth [13], answering a question of von Neumann; see also 
Pippenger [18]). Here we see a similar phenomenon when there are no carries. In studying 
Conjectures 3 and 4, we observed that the numbers of divisors for the runners-up, 2 k — 3 and 
2 k — 2 k ~ 2 — 1, appeared to be converging to one-fifth of the number of divisors of 1 1 ... 1 1 1 2 • 
This is proved in Theorem 19. Our proof is modeled on Knopfmacher and Robbins's proof 
[12] of (23), and we present the proof in such a way that it yields both results simultaneously. 

The sum-of-di visors function a^n) is the subject of §7. There are analogues of the 
perfect numbers, although they seem not to be as interesting as in the classical case. Section 
8 discusses the dismal analogue of the partition function. Theorems 22 and 23 give explicit 
formulas for the number of partitions of n into distinct parts. 

This is the second of a series of articles dealing with various kinds of carryless arithmetic, 
and contains a report of our investigations into dismal arithmetic carried out during the 
period 2000-2011. This work had its origin in a study by the second author into the results 
of performing binary arithmetic calculations with the usual addition and multiplication of 



binary digits replaced by other operations. If addition and multiplication are replaced by the 
logical operations OR and AND, respectively, we get base 2 dismal arithmetic. (If instead 
we use XOR and AND, the results are very different, the squares for example now forming 
the Moser-de Bruijn sequence A000695.) Generalizing from base 2 to base 10 and then to 
an arbitrary base led to the present work. 

In the first article in the series, [1], addition and multiplication were carried out "mod 
10", with no carries. A planned third part will discuss even more exotic arithmetics. 

Although dismal arithmetic superficially resembles "tropical mathematics" [20], where 
addition and multiplication are defined by x © y := min{x,y}, x y := x + y, there is no 
real connection, since tropical mathematics is defined over 1U {oo}, uses carries, and is not 
base-dependent. 

Notation. The base will be denoted by b and the largest digit in a base b expansion by 
/3 := b — 1. We write n = nk-iTik-2 ■ • •fta7io|& to denote the base b representation of the 
number Yli=o n $ \ anc ^ we define len&(n) := k. The components rij will be called the digits 
of n, even if b ^ 10. In the examples in this paper b will be at most 10, so the notation 
rik-i ■ ■ ■ riinolb (without commas) is unambiguous. The symbols +b and x b denote dismal 
addition and multiplication, and we omit the base b if it is clear from the context. All 
non-bold operators (+, x, <, etc.) refer to ordinary arithmetic operations, as do unqualified 
terms like "smallest," "largest," etc. We usually omit ordinary multiplication signs, but 
never dismal multiplication signs. We say that p divides n in base b (written p -<& n) if 
px b q = n for some q, and that p = Pk-i ■ ■ ■ PiPo\b is dominated by n = rik-i ■ ■ ■ nirio\b 
(written p <C& n) if pi < rij for all %. The symbol "];," always marks the end of a base b 
expansion of a number, and is never used for "divides in base 6." 

2 Basic definitions and properties 

We began, as we all did, in base 10, but from now on we will allow the base b to be an 
arbitrary integer > 2. 

Let A denote the set of base b "digits" {0, 1, 2, . . . , b — 1}, equipped with the two binary 
operations 

m+t,n := max{m, n}, mX b n := min{m, n}, ioim,n^A. (5) 

A dismal number is an element of the semiring A[X] of polynomials J^i=o n %^-' \ n i ^ A. If 
M[X] := ^;=o rn i^- t an d N[X\ := X]j=o n *^ 4 are dismal numbers then their dismal sum is 
formed by taking the dismal sum of corresponding pairs of digits, analogously to ordinary 
addition of polynomials: 

max{fc,/}— 1 

M[X}+ b N[X] := J2 P* Xl > ( 6 ) 

j=0 

where pi := m, +b n ii and their dismal product is similarly formed by using dismal arithmetic 
to convolve the digits, analogously to ordinary multiplication of polynomials: 

k+l-2 

M[X]x b N[X] := Yl ^' (7) 

j=0 
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We will identify a dismal number N(X) = Yli=o n iX % with the integer n whose base b 
expansion is n = Yli=o n i^ ( n * s obtained by evaluating the polynomial N(X) at X = b), 
and we define len b (n) := k. The rules (5)-(7) then translate into the rules for dismal addition 
and multiplication stated in §1: there are no carries, and digits are combined according to 
the rules in (5). The fe-ary number ^i=o n & % wm a ^ so De written as Uk-iUk-2 ■ ■ ■ n\no\ b . 
Dismal numbers are, by definition, always identified with nonnegative integers. Note that 
len b (m+ b n) = max{lenft(m),lenb(n)} and len&(mX&n) = len&(m) + len&(n) — 1. 

Theorem 1. The dismal operations + b and x b on A[X] satisfy the commutative and asso- 
ciative laws, and x b distributes over + b . 

Proof. (Sketch.) Each law requires us to show the identity of two polynomials, and so reduces 
to showing that certain identities hold for the coefficients of each individual degree in the 
two polynomials. These identities are assertions about min and max in the set A, which hold 
since (A, <) is a totally ordered set, and (A, min, max) is a distributive lattice (cf. [9]). □ 

If R denotes the operation of reversing the order of digits and m and n have the same 
length, then R{m + b n) = R{m) + b R(n) and R{m X b n) = R{m) X b R{n). 

Individual digits in a dismal sum or product can often be varied without affecting the re- 
sult, so dismal subtraction and division will not be defined. Example: 16 +i 75 = 26 +io 75 = 
76, 16 x 10 75 = 16 x 10 85 = 165. This is why dismal numbers form only a semiring. On the 
other hand, this semiring does have a multiplicative identity (see the next section), and there 
are no zero divisors. 

In certain situations we can give a more precise statement about how changing digits in 
a dismal sum or product affects the result. We begin with a lemma about ordinary functions 
of real variables. 

Lemma 2. Let f be a single-valued function of real variables X\, . . . ,Xk, k > 2, formed by 
repeatedly composing the functions (x,y) !->■ mm{x,y} and (x,y) i-> max{x,y}. If g is a 
nondecreasing function of x, meaning that 

x < y => g(x) < g(y) , (8) 

then 

f(g(xi),...,g(x k )) = g(f(x 1 ,...,x k )) 1 (9) 

for all real xi, . . . ,xj,- 

We omit the easy inductive proof. 

We define a base b digit map to be a nondecreasing function g mapping {0, 1, . . . , b—1} into 
itself. The map g need not be one-to-one or onto. If g is a digit map and n = n^-i . . . nino\ b 
then we set g(n) := g(n k -i) ■ . ■g(n 1 )g(n )\ b . 



Theorem 3. If m and n are dismal numbers and g is a base b digit map, then 

g(m+ b n) = g(m)+ b g(n), 

g(mx b n) = g(m)x b g{n). (10) 

Proof. This follows from Lemma 2, since the individual digits of m+ b n and mx b n are 
functions of the digits of m and n of the type considered in that lemma. □ 

Corollary 4. Ifp = mx b n then p can also be written as m! X b n' ' , where m' and n' use only 
digits that are digits of p. 

Proof. (Sketch.) Arrange all distinct digits occurring in p, m, and n in increasing order. 
Then construct a digit map g by increasing or decreasing the digits in m and n that are not 
in p until they coincide with digits of p, leaving the digits of p fixed. □ 

For example, consider the product 165 = 16 x 10 85 mentioned above. The digits involved 
are 1, 5, 6, 8, and the digit map described in the proof fixes 1, 5, and 6 and maps 8 to 6. The 
resulting factorization is 165 = 16 x 10 65. (The additive analogue of Corollary 4 is true, but 
trivial.) 

We will see other applications of Theorem 3 in the next section. 

Note that when we are computing the base b dismal sum or product of two numbers p 
and q, once we have expressed p and q in base b, the value of b plays no further role in the 
calculation. Of course we need to know b when we convert the result back to an integer, but 
otherwise b is not used. So we have: 

Lemma 5. If the largest digit that is mentioned in a base b dismal sum or product is d 
(where < d < b — 1), then the same calculation is valid in any base that exceeds d. 

For example, here is the calculation of the base 2 dismal product of 13 = 1101 12 and 
5 = 101| 2 : 
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This tells us that 13 x 2 5 = 61, but the same tableau can be read in base 3, giving 37 x 3 10 = 
361, or in base 10, giving 1101 x 10 101 = 111101. 

Recall that we say that p divides n in base b (written p -< b n) if p x b q = n for some q. 
Since \en b (p) < len b (n), nonzero numbers have only finitely many divisors. We also say that 
p := pk-i ■ ■ -PiPo\b is dominated by n := rik-i . . . nino\ b (written p <Cb n) if pi < n t for all i. 
Then p <C;, n if and only if p + b n = n. Another consequence of Lemma 2 is: 

Lemma 6. If p <C& m and q <t^ b n then p+ b q <ti b m+ b n and px b q <^ b rn x b n. 

Finally, we remark without giving any details that, in any base, the sets of numbers with 
digits in nondecreasing order, or in nonincreasing order (see A009994 and A009996 for base 
10) are closed under dismal addition and multiplication. 
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3 Dismal primes 

In dismal arithmetic in base b, for bases b > 2, the multiplicative identity is no longer 1 (for 
example, 1 x 10 23 = 11, not 23). In fact, it follows from the definition of multiplication that 
the multiplicative identity is the largest single-digit base b number, (3 := b — 1. For base 
b = 10 we have = 9, and indeed the reader will easily check that 9x 10 n = n for all n. An 
empty dismal product is defined to be j3, by convention. 

If p Xf, q = (3, then p = q = (3, so (3 is the only unit. We therefore define a prime in base b 
dismal arithmetic to be a number, different from (3, whose only factorization is (3 times itself. 

If p is prime, then at least one digit of p must equal [3 (for if the largest digit were r < (3, 
then p = r x b p, and r would be a divisor of p). The base b expansions of the first few primes 
are 1(3 (this is the smallest prime), 2/3, 3/3 ... (3 — 1 (3, (30, (31, . . . , (3/3, 10/3, .... In base 10, the 
primes are 

19, 29, 39, 49, 59, 69, 79, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 109, 209, 219, 

309, 319, 329, 409, 419, 429, 439, 509, 519, 529, 539, 549, 609, 619, 629, 639, . . . (11) 

(A087097). Notice that the presence of a digit equal to (3 is a necessary but not sufficient 
condition for a number to be a prime: 11/5 1& = 1/3 1& x b l/3|& is not prime (see A087984 for 
these exceptions in the case b = 10). In base 2, the primes (written in base 2) are 

10, 11, 101, 1001, 1011, 1101, 10001, 10011, 10111, 11001, 11101, 100001, ... (12) 

(A171000). In view of the interpretation of base 2 dismal arithmetic in terms of Boolean 
operations mentioned in §1, the corresponding polynomials 

X, X + 1, X 2 + 1, X 3 + 1, X 3 + X + 1, X 3 + X 2 + 1, A 4 + 1, X 4 + X + 1, . . . , 

together with 1, might be called the OR- irreducible Boolean polynomials. Their decimal 
equivalents, 

1, 2, 3, 5, 9, 11, 13, 17, 19, 23, 25, 29, 33, 35, 37, . . . 

form sequence A067139 in [17], contributed by Jens VoB in 2002. 

All numbers of the form 100 . . . 00/3 \b (with zero or more internal zeros) are base b primes, 
since there is no way that pX b q can have the form 100 . . . 00/3 1& unless p or q is a single-digit 
number. So there are certainly infinitely many primes in any base. 

Since 1/3 \b is the smallest prime, the numbers 1/5 1 1, x b n are another analogue of the even 
numbers. Whereas the first version of the even numbers, given in (1) for base 10, was simply 
"replace all digits in n that are bigger than 2 with 2's," this version is more interesting. In 
base 10 we get 

0, 11, 12, 13, 14, 15, 16, 17, 18, 19, 110, 111, 112, 113, 114, 115, 116, 117, 118, . . . (13) 

(A162672, which contains repetitions and is not monotonic). 

If b > 2, there are numbers which cannot be written as a product of primes (e.g., 1). 

Theorem 7. Any base b number with a digit equal to (3 is a (possibly empty) dismal product 
of dismal primes. 



Proof. The number (3 itself is the empty product of primes. Every two-digit number with /3 
as a digit is already a prime. If there are more than two digits, either the number is a prime, 
or it factorizes into the product of two numbers, both of which must have /3 as a digit. The 
result follows by induction. □ 

Corollary 8. Every base b number can be written as r times a dismal product of dismal 
primes, for some r e {0, 1, 2, . . . , b — 1}. 

Proof. Let r be the largest digit of n. If r = /3 the result follows from the theorem. Otherwise, 
let m be obtained by changing all occurrences of r in the 6-ary expansion of n to /3's, so that 
n = r x&m, and apply the theorem to m. □ 

Even when it exists, the factorization into a dismal product of dismal primes is in gen- 
eral not unique. In base 10, for example, the list of numbers with at least two different 
factorizations into a product of primes is 

1119, 1129, 1139, 1149, 1159, 1169, 1179, 1189, 1191, 1192, 1193, 1194, 1195, . . . , (14) 

where for instance 1119 = 19 x 10 19 x 10 19 = 19 x 10 109 (A171004). 
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Table 1: Numbers of dismal primes with k digits in bases 2 and 10 (A169912, A087636). 

We can, of course, study the primes dividing n, even if n does not contain a digit equal 
to (3. Without giving any details, we mention that [17] contains the following sequences: the 
number of distinct prime divisors of n (A088469), their dismal sum (A088470), and dismal 
product (A088471); 2 also the lists of numbers n such that the dismal sum of the distinct 
prime divisors of n is < n (A088472), < n (A088473), > n (A088475), > n (A088476); as 
well as the numbers n such that the dismal product of the distinct prime divisors of n is < n 
(A088477), < n (A088478), = n (A088479), > n (A088480), and > n (A088481). There is 
no analogue of A088476 or A088481 in ordinary arithmetic. 

One omission from the above list is explained by the following theorem. 



2 A088471 has an unusual beginning: 9, 9, 9, 9, 9, 9, 9, 9, 9, 90, 123456789987654321, 19, 19, 19, . 



Theorem 9. In base b dismal arithmetic, n is prime if and only if the dismal sum of its 
distinct dismal prime divisors is equal to n. 

Proof. If n = p is prime then the sum of the primes dividing it is p. Suppose n is not prime 
and let m be the sum of the distinct dismal primes diving n. If n is divisible by a prime 
p with \en b (p) = len b (n), then n = r x b p, r < (3, the largest digit in n is r, and som/n 
since n <C& p <Cft m. If len.b(p) < len&(n) for all prime divisors p, then len^ra) < len&(n), and 
again m ^ n. D 

We now consider how many primes there are. Let iTb(k) denote the number of base b 
dismal primes with k digits. Table 1 shows the initial values of ^(k) and 7Tio(/c). Necessary 
conditions for a number n to be prime are that it contain (3 as a digit and (if k > 2) does 
not end with 0. There are 

(b - l) 2 b k ~ 2 - (b - 2) (b - if- 2 (15) 

such numbers. It seems likely that, as k increases, almost all of these numbers will be prime, 
and the data in Table 1 is consistent with this. We therefore make the following conjecture. 

Conjecture 1. 

n k (b) ~ (b-l) 2 b k ~ 2 ask^oo. (16) 

We can get a lower bound on Hb{k) by producing large numbers of primes, using the 
process of "promotion." We call a base b number with at least two digits a pseudoprime if 
its only factorizations are of the form n = pXbl where at least one of p and q has length 
1. In base b, n is a prime if and only if it is a pseudoprime and contains a digit 0. If 
nfc_infc_2...no|& is a pseudoprime and r is its maximal digit, then nfc_irifc„ 2 ---^o|r+i is a base 
r + 1 prime and furthermore nk-\nk-2---no\ c is a pseudoprime in any base c > r + 1. In base 
2 there is no difference between primes and pseudoprimes. As long as we exclude numbers 
ending with 0, reversing the digits of a number does not change its status as a prime or 
pseudoprime. 

The advantage of working with pseudoprimes rather than primes is that the inverse image 
of a pseudoprime under a digit map (see §2) is again a pseudoprime. 

Theorem 10. For a base b digit map g, if g{nk_ink-2---no\b) is a pseudoprime and g(rik-i) 

is not 0, then nk-\nk-2---no\b is a pseudoprime. 

Proof. This follows immediately from Theorem 3. □ 

So if p is a pseudoprime, then any number n with the property that there is a digit map 
sending n to p is also a pseudoprime; we think of n as being obtained by "promoting" p, and 
call p the "template" for n. 

Here is an equivalent way to describe the promotion process. Suppose the distinct digits 
in p, the template, or number to be promoted, are d\ < d 2 < d 3 , . . .. For each d iy choose a 
set of digits S(di) such that all the digits in S(di) are strictly less than all those in S(dj), 
for % < j. Replace any digit di in p by any digit in S(di). All numbers n obtained in this 
way are promoted versions of p (the required digit map g being defined by g(c) = di for all 

ceS{di)). 

9 



Any pseudoprime (in any base) with at most four digits can be obtained by promoting 
a base 2 prime. At length 2, 11| 2 is a prime (see (12)), so every 2-digit number rs\b is a 
pseudoprime, using the digit map g that sends r and s to 1. This is valid even if s = 0, 
since (8) still holds. At length 3, 101 1 2 is a base 2 prime, so any three-digit number rst\b 
with r > s, t > s is a pseudoprime (take g{r) = g(t) = l,g(s) = 0), and this captures all 
three-digit pseudoprimes. There are three templates of length 4, 1001|2, 1011 1 2 , and 1101 1 2- 
These can be promoted to capture all four-digit primes, which are the numbers rstu\b for 
which one of the following holds: 

r and u are both strictly greater than s and t , 
r, s, and u are strictly greater than t , 
r, t, and u are strictly greater than s . 

For lengths greater than four, we must use some nonbinary templates to capture all 
pseudoprimes, and as the length k increases so does the fraction of nonbinary templates 
required, as shown in Table 2. The columns labeled (a) and (b) give the number of binary 
templates and the total number of templates, respectively, and the columns (c) and (d) give 
the number of base 10 primes obtained by promoting the templates in columns (a) and (b). 



k 


(a) 


(b) 


(c) 


(d) 


2 


1 


1 


18 


18 


3 


1 


1 


81 


81 


4 


3 


3 


1539 


1539 


5 


5 


8 


17661 


20457 


6 


9 


51 


135489 


242217 



Table 2: For lengths k = 2 through 6, the numbers of (a) binary templates, (b) all templates, 
(c) base 10 primes obtained by promoting the binary templates, and (d) base 10 primes 
obtained by promoting all the templates. 

We can reduce the list of templates by omitting those that are reversals of others. Table 
3 shows the reduced list of templates of lengths < 6. 

Since 1 00 . . . 1| 2 is prime, the promotion process tells us for example that the numbers 

fe-2 

r S1S2 ■ ■ ■ s/c-2 P\b are prime, for 1 < r < (3, provided each of the digits Si is in the range 

fc-2 

through r — 1. This gives (b — l) k ~ 2 + 2(6 — 2) k ~ 2 -\ — • = 0(b k ~ 2 ) primes of length k, which 
for b > 2 is of exponential growth but smaller than (16). 

4 Dismal squares 

One might expect that it would be easier to find the number of dismal squares of a given 
length than the number of dismal primes, but we have not investigated squares as thoroughly, 



10 



11 


100001 


102212 


120212 


101 


100011 


102221 


120221 


1001 


100101 


103223 


120222 


1011 


100111 


103233 


121022 


10001 


101011 


110212 


121102 


10011 


101221 


112021 


122102 


10111 


101222 


112022 


122202 


12021 


102201 


120021 


132023 


12022 


102202 


120022 


133023 



Table 3: Reduced list of prime templates: every pseudoprime of length < 6 can be obtained 
by promoting one of these 36 pseudoprimes or its reversal (A191420). 

and we do not even have a precise conjecture about their asymptotic behavior. In base 2, 
the first few dismal squares, written in base 10, are 



0, 1, 4, 7, 16, 21, 28, 31, 64, 73, 84, 95, 112, 125, 124, 127, 256, 273, . . . 



;i7) 



(A067398, also contributed by Jens VoB in 2002), and the numbers of squares of lengths 1 
(including 0), 3, 5, 7, . . . are 

2, 2, 4, 8, 15, 29, 55, 105, 197, 367, 678, 1261, 2326, 4293, 7902, 14431, . . . (18) 

(A190820). In base 10, the first few squares were given in (2), and the numbers of squares 
of lengths 1 (including 0), 3, 5, 7, . . . are 

10, 90, 900, 9000, 74667, 608673, . . . 

(A172199). The sequence of base 10 squares is not monotonic (for example 1011 < 1020 
yet 1011 x 10 1011 = 1011111 > 1020 x 10 1020 = 1010200), and contains repetitions. The 
numbers which are squares in more than one way are 

111111111, 111111112, 111111113, 111111114, 111111115, 111111116, 111111117, ... , 

e.g., 111111111 = 11011 x 10 11011 = 11111 x 10 11111 (A180513, A181319). 

We briefly mention two other questions about squares to which we do not know the 
answer: (i) In base 2, how many square roots does 2 2k+1 — 1 have? This is a kind of 
combinatorial covering problem. For k — 0, 1, . . . the counts are 

1, 1, 1, 1, 2, 3, 5, 9, 15, 28, 50, 95, 174, 337, 637, 1231, 2373, 4618, 8974, 17567, 34387, 
67561, 132945, 262096, 517373, 1023366, 2025627, 4014861, 7964971, 15814414, 
31424805, 62490481, 124330234, 247514283, 492990898, 982307460, 1958093809, 
3904594162, 7788271542, 15539347702, 31012331211, . . . (19) 

(A191701). Is there a formula or recurrence for this sequence? (ii) In base b, if we consider 
all p such that px b p = n, does one of them dominate all the others (in the 3>& sense)? If so 
the dominating one could be called the "principal" square root. 
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5 The divisibility poset 

One drawback to dismal arithmetic is that there is more than one way to order the dismal 
numbers, and no ordering is fully satisfactory. 

The usual order on the nonnegative integers (< or <) is unsatisfactory, since (working 
in base 10) we have 18 < 25 yet 18 + 32 = 38 > 25 + 32 = 35, 32 < 41 yet 32 x 3 = 32 > 
41x3 = 31. 

The dominance order (<Cb) is more satisfactory, in view of Lemma 6 and the distributive 
law of Theorem 1, but has the drawback that m divides n does not imply that m <Cb n (e.g., 
12|i divides 11 |io, yet lljio <io 12| 10 ). 

The partial order induced by divisibility (~<b) is worth discussing, as it has some interest- 
ing properties and is the best way to look at dismal numbers as long as we are considering 
only questions of factorization and divisibility. For simplicity we will restrict the discussion 
to base 10. 



11 
10 

9 

8 
7 

4 
3 
2 
1 





Figure 1: Beginning of the divisibility poset (see text for description). The left-hand column 
gives the rank (A161813). 

Figure 1 displays the beginning of the Hasse diagram ([21, p. 99]) of this partially ordered 
set (or poset), and shows all positive numbers with one or two digits, and a few larger 
numbers. There are too many edges to draw in the diagram, so we will describe them in 
words. 
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The multiplicative identity 9, the "zero element" in the poset, is at the base. The other 
single-digit numbers 8-<7-<6-<----<l are above it in the left-hand column, and 8 is 
joined to 9. The numbers are arranged in rows according to their rank (shown at the extreme 
left of the diagram). The numbers of rank 1 consist of 8 and the (infinitely many) primes: 
19, 29, ... , 91, 90, 109, 209, 219, 309, . . . (A144171). All of these are joined to 9. The numbers 
of rank 2 are 7, 18, 28, ... , 81, 80, 108, 119, . . . (A144175), and so on. Every two-digit number 
of rank h (1 < h < 9) is joined to the single-digit number h — 1 on the left of the diagram., 
as indicated by the square brackets. A two-digit number to the left of the central column of 
the pyramid is joined to the number diagonally below it to the left (e.g., 56 is joined to 57), 
and a two-digit number to the right of the central column is joined to the number diagonally 
below it to the right (e.g., 65 is joined to 75). A number in the central column is joined 
to the three numbers immediately below it (e.g., 66 is joined to 67,77,76). A number rO 
(1 < r < 9) in the right-hand column is joined to the number immediately below it and to 
the single- digit number r. 

Only a few numbers with more than two digits are shown, but a more complete diagram 
would show for example that 1 is joined to, besides 10 and 11, many other decimal numbers 
whose digits are 0's and l's, such as 101, 1001, 1011, . . .. all of rank 9. The figure is complete 
in the sense that all downward joins are shown for all the numbers in the diagram. 

One perhaps surprising property of the divisibility poset is that the greatest lower bound 
(or greatest common divisor) m An and the least upper bound (or least common multiple) 
mVnof two dismal numbers m and n need not exist, and so this poset fails to be a lattice 
[9, Chap. 1]. For example, again working in base 10, the rank 2 numbers 8989 and 9898 
are each divisible by (and joined to) the nine primes 909, 919, . . . , 989. However, these nine 
primes are incomparable in the -<io order, so neither 8989 A 9898 nor 909 V 919 exist. 

Although greatest common divisors need not exist, we can still define two dismal numbers 
to be relatively prime if their only common divisor is the unit j3. 

In the next section we will study the number of divisors function db(n). Candidates for 
divisors of n are all numbers m with len^m) < lenfe(n). It is therefore appropriate to define 
the dismal analogue of the Euler totient function, 0&(n), to be the number of numbers m 
with lenb(m) < len;,(n) which are relatively prime to n. The initial values of 02 {n) and 4>io(n) 
are shown in Table 4. 



n 


12 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 


<h(n) 


1 2 2 4 6 2 4 8 14 6 14 5 14 5 7 16 30 14 30 12 


01OW 


11111111 9 18 2 18 18 18 18 18 18 18 90 18 



Table 4: Values of totient functions fain) and 0io(^) (A191674, A191675). 



6 The number of dismal divisors 

Let db(n) denote the number of dismal divisors of n in base b, and let <Jb(n) denote the 
dismal sum of the dismal divisors of n. These functions are more irregular than their classical 
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analogues, as can be seen from the examples in Table 5, and there are no simple formulas for 
them. In this section we study some of the properties of d&(n). Note that if r is the largest 
digit in n, then the smallest divisor of n is r, and the largest divisor is the number obtained 
by changing all the r's in n to /3's. 



n 


divisors (base 10) 


d w (n) 


o-io (n) 


1 


1,2,3,4,5,6,7,8,9 


9 


9 


2 


2,3,4,5,6,7,8,9 


8 


9 


3 


3,4,5,6,7,8,9 


7 


9 


4 


4,5,6,7,8,9 


6 


9 


5 


5,6,7,8,9 


5 


9 


6 


6,7,8,9 


4 


9 


7 


7,8,9 


3 


9 


8 


8,9 


2 


9 


9 


9 


1 


9 


10 


1,...,9, 10,20, ...,90 


18 


99 


11 


1, . . . , 9, rs with 1 < r, s < 9 


90 


99 


12 


2,..., 9, 12, 13,..., 19 


16 


19 



Table 5: In base 10, the dismal divisors of the numbers 1 through 12 and the corresponding 
values of d w (n) and a w (n) (A189506, A087029, A087416). 

A base b dismal prime p has two divisors, 6—1 and p, so db(p) = 2. In the other direction, 
a divisor of a fc-digit number n has at most k digits, so 



2 < d b (n) < b k - 1 . 



(20) 



Base b numbers of the form 111 . . . 1|& (that is, in which all the base b digits are 1) come 
close to meeting this upper bound — see Remark (iv) following Theorem 13. We make the 
following conjectures. 

Conjecture 2. In any base b > 3, among all k-digit numbers n, db(n) has a unique maximum 
atn = (b k -l)/(b-l) = 111 . . . lj 6 . 

Conjecture 3. In base 2, among all k-digit numbers n, the maximal value of c?2(n) occurs 
at n = 2 k — 2 = 111 . . . 1 1 2 , and this is the unique maximum for n ^ 2,4. 

Conjecture 4. In base 2, among all odd k-digit numbers n, d2(n) has a unique maximum 
at n = 2 k — 1 = 111 ... Ill 1 2, and if k > 3 and k 7^ 5, the second- largest value ofd 2 (n) occurs 
at n = 2 k - ; 
n. 



111... 101k n = 2 fe -2 



fc-2 



101 ... 111k an d possibly other values of 



The numerical evidence supporting these conjectures is compelling. For example, in base 
10, if we study the sequence d w (n), n > 1 (A087029) for n < 10 6 , and write down di (n) 
each time it exceeds dio(m) for all m < n, we obtain the values 

9, 18, 90, 180, 819, 1638, 7461, 14922, 67968 , 



11 



(see A186443) at these (decimal) values of n: 

1, 10, 11, 110, 111, 1110, 1111, 11110, 11111 . 

If n is a 5-digit decimal number, the eight largest values of dio{n) are, in decreasing order, 

67968, 39624, 21812, 14922, 11202, 9616, 6732, 6570, 

at these values of n: 

11111, 22222, 33333, 11110,44444, 12222 or 22220 or 22221, 11011, 10111 or 11101. 

The number n = (10 fc — 1)/9 = 111 . . . 1| 10 is a clear winner among all fc-digit decimal numbers 
for k < 5. The data for bases 3 through 9 is equally supportive of Conjecture 2. Likewise, the 
binary data strongly supports Conjectures 3 and 4 — see Table 6 for the initial values of d 2 {n). 
Table 6 also shows why k = 5 is mentioned as an exception in Conjecture 4: among 5-digit 
odd numbers, ^(llOll^) = 4 is the runner-up, ahead of o?2(10111|2) = ^(HlOl^) = 2. 



n (base 2) 


d 2 {n) 


n (base 2) 


d 2 (n) 


n (base 2) 


d 2 (n) 


n (base 2) 


d 2 (n) 


— 


— 


1000 


4 


10000 


5 


11000 


8 


1 


1 


1001 


2 


10001 


2 


11001 


2 


10 


2 


1010 


4 


10010 


4 


11010 


4 


11 


2 


1011 


2 


10011 


2 


11011 


4 


100 


3 


1100 


6 


10100 


6 


11100 


9 


101 


2 


1101 


2 


10101 


3 


11101 


2 


110 


4 


1110 


6 


10110 


4 


11110 


10 


111 


3 


mi 


5 


10111 


2 


11111 


8 



Table 6: In base 2, the number of dismal divisors of the numbers 1 through 31 (A067399). 

For a more dramatic illustration of Conjectures 2 and 3, see the graphs of sequences 
A087029 and A067399 in [17]. Although the evidence is convincing, we have not, unfortu- 
nately, succeeded in proving these conjectures. 

We are able to determine the exact values of d b (lll . . . Ill |^) for all b (the conjectural 
winner for b > 3 and the conjectural winner among odd numbers in the binary case) and 
d 2 (lll . . . 101 12) = g?2(101 . . . Ill) (the conjectural runners-up among odd binary numbers of 
length k > 5). 

We begin with a lemma that describes the effect of trailing zeros. 

Lemma 11. // the base b expansion of n ends with exactly r > zeros, so that n = mb r , 
with b \ m, then 

d b (n) = (r + V)d b {m) . 

Proof, li px b q = m then b j p, b j q, and the r + 1 numbers pb % (0 < i < r) dismally divide 
n, since 

{pV)x b {qb r - l )=n. 

Conversely, if p' x b q' = n then p' = pb 1 , q' = qb r ~ l , b \ p, b \ q, for some i with < i < r. So 
each dismal divisor of m corresponds to exactly r + 1 dismal divisors of n. □ 
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For example, in base 10 the dismal divisors of 7 are 7, 8, 9 and the dismal divisors of 700 
are 7, 8, 9, 70, 80, 90, 700, 800, 900. 

One reason Conjectures 2-4 seem hard to prove is the erratic behavior of db(n). In contrast 
to the above lemma, the effect of internal zeros is hard to analyze. Suppose the i-th digit in 
the fc-ary expansion of n is zero. This implies that if n = p x&g, then all entries in the i-th 
column of the long multiplication tableau must be zero, which imposes many constraints on 
the 6-ary expansions of p and q. One would expect, therefore, that changing the zero digit 
to a larger number — thus weakening the constraints — would always increase the number of 
divisors of n. Roughly speaking, this is true, but there are many cases where it fails. For 
example, <i 2 (lllll| 2 ) = 8, but <i 2 (11110| 2 ) = 10 (see Table 6, Lemma 11 and Theorem 12). 
Again, d 2 (10101| 2 ) = 3, but 10111 1 2 is prime, so ci 2 (10111| 2 ) = 2. In any base b, 11 00 ... |& 

r 

has (r + l)((b — l) 2 + (b — 1)) divisors, whereas 11 00 . . .0 1|& has (b — l) 3 + (6 — 1) divisors, 

r-1 

a smaller number if r is large. 

The next result was conjectured by the second author and proved by Richard Schroeppel 
in 2001 [15]. An alternative proof (via a bijection with a certain class of polyominoes) was 
given by Frosini and Rinaldi in 2006 [6]. We give a version of Schroeppel's elegant direct 
proof, partly because it has never been published, and partly because we will use similar 
arguments later. 

Theorem 12. In base 2, the number of dismal divisors of 111 . . . 1 j 2 (with k 1 's) is equal to 
the number of compositions of k into parts of which the first is at least as great as all the 
other parts. 

Proof. Suppose p*iq = lll...l|2 (with k l's) where len 2 (p) = r, len 2 (g) = k + 1 — r. 
By examining the long multiplication tableau for px 2 q, we see that it is also true that 
P x 2 q' = HI • • • 1 1 2 where q' = 111 . . . 1| 2 , with k + l—r l's (for if there is a 1 in each column 
of the tableau for px 2 g, that is still true for px 2 g'). So in order to find all the divisors p 
of 111 ... 1 12 (with k l's) we may assume that the cofactor q has the form 111 . . . 1| 2 (with s 
l's, for some s, 1 < s < k). 

We establish the desired result by exhibiting a bijection between the two sets. Let 
k = c\ + c 2 + . . . + ct be a composition of k in which c\ > q > 1 (2 < % < t). Let ip(ci) denote 
the binary vector 000 ... 01 with q — 1 0's and a single 1. The divisor p corresponding to 
this composition has binary representation given by the concatenation 

1^(C 2 )^(C 3 ) ...lf>(Cr)\ 2 , (21) 

of length k + 1 — C\. If we set q = 111 . . . 1| 2 , of length c±, then p X 2 q = 111 . . . 1| 2 (with 
k l's). This follows from the fact that if the binary representation of p contains a string of 
exactly s 0's: 

... 1 000 ... 1 ... , 

s 

then, when we form the product p X 2 g, the 1 immediately to the right of these 0's will 
propagate leftward to cover the 0's if and only if len 2 (g) > s + 1, which is exactly the 
condition that C\ > Ci for all i > 2. □ 
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For example, the eight compositions of 5 in which no part exceeds the first and the 
corresponding factorizations p X 2 q of 11111J2 are shown in Table 7. Dots have been inserted 
in p to indicate the division into the pieces ip{cj). 



Composition of 5 


divisor p 


cofactor q 


5 


1| 2 


11111|2 


41 


1.1| 2 


1111| 2 


32 


i.oi| 2 


Hl| 2 


311 


1.1. 1| 2 


Hl| 2 


221 


1.01. 1| 2 


H| 2 


212 


1.1.01| 2 


H| 2 


2111 


1.1.1.1I2 


H| 2 


11111 


1.1.1.1.1I2 


1| 2 



Table 7: Illustrating the bijection used to prove Theorem 12. 

Remarks. 

(i) Using the bijection defined by (21), the number of l's in the binary expansion of the 
divisor p is equal to the number of parts in the corresponding composition. 

(ii) It follows immediately from the interpretation in terms of compositions that the 
numbers d^i^ — 1) have generating function 



00 

k=\ 


i(2 k - 


-l)z k -- 


00 
1=1 


-(* 


z l 
+ z 2 + - 


■ + z l ) 








00 
- V 


(1- 


■ z)z l 







2z + z 



1+1 



(22) 



(the index of summation, I, corresponds to the first part in the composition). 

(iii) The initial values of this sequence are shown in Table 8. This sequence appears in 



k 


1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 


d 2 (2 k - 1) 


1 2 3 5 8 14 24 43 77 140 256 472 874 1628 3045 5719 



Table 8: Values of d 2 (2 k - 1). 

entries A007059 and A079500 in [17], although the indexing is different in each case. The 
sequence also occurs in at least four other contexts besides the two mentioned in Theorem 
12, namely in the enumeration of balanced ordered trees (Kemp [10]), of polyominoes that 
tile the plane by translation (Beauquier and Nivat [2], Brlek et al. [3]), of Dyck paths (see 
A007059), and in counting solutions to the postage stamp problem (again see A007059). 
The article by Frosini and Rinaldi [6] gives bijections between four of these six enumerations. 
In this context we should also mention the recent article of Rawlings and Tiefenbruck [19], 
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which, although not directly related to the problems we consider, discusses other connections 
between the enumeration of compositions, permutations, polyominoes, and binary words. 

(iv) The asymptotic behavior of this sequence is quite subtle. From the work of Kemp 
[10] and Knopfmacher and Robbins [12] it follows that 



d 2 (2 h 



Hog 2 



'1 + Gfc), as k — > oo , 



(23) 



where 6^ is a bounded oscillating function with |Ofc| < 10 5 (see the proof of Theorem 19 

below). 

In order to determine d&(lll . . . 1|&) for bases b > 2, we first classify compositions in 
which no part exceeds the first according to the number of parts. Let T(k,t) denote the 
number of compositions of k into exactly t parts (with 1 < t < k) such that no part exceeds 
the first. Table 9 shows the initial values. This is entry A184957 in [17]. 3 



k\t 


12 3 4 


5 6 7 8 


1 


1 




2 


1 1 




3 


1 1 1 




4 


12 11 




5 


12 3 1 


1 


6 


13 4 4 


1 1 


7 


13 6 7 


5 1 1 


8 


1 4 8 11 


11 6 1 1 



Table 9: Initial values of of T(k, t), the number of compositions of k into exactly t parts such 
that no part exceeds the first. 

The values of T(k, t) are easily computed via the auxiliary variables ^(k, t, m), which we 
define to be the number of compositions of k into t parts of which the first part, m, is the 
greatest (for 1 <t < k, 1 < m < k). We have the recurrence 



min{m,fc+2— t— m} 

nf(k,t,m) = Yl l(k-j,t-l,m) 



(24) 



(classifying compositions according to the last part, j), for m > 1, t > 1, t + m < k — 1, with 
initial conditions 

j(k,t,l) = 6 t ,k, 
~f(h,l,m) = 5 mtk , 



3 An array equivalent to this, A156041, was contributed to [17] by J. Grahl in 2009 and later studied by 
A. P. Heinz and R. H. Hardin. 
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where d~ij = 1 if % = j or if % 7^ j . Then 

k+l-t 

T(k,t) = J2 l(k,t,m). (25) 

m=l 

Since j(k, t, m) is the coefficient of z k in z m (z + z 2 + ■ • • + z m ) t_1 , it follows that column t 
of Table 9 has generating function 

EnM) Z ^ ir ^EMr( i ; 1 rf^r- (*» 

fe=l V y r=0 V 7 

Since the total number of compositions of k into t parts is ( t Zi)> an d in at least a fraction 
I of them the first part is the greatest, we have the bounds 

Theorem 13. 

db {j~r) = MU_^J,\b) = J2T(k,t)(b-lf. (28) 

k t=l 

Proof. Suppose p x b q = 11 ... 1 1 6 . At least one of p and q, say q, must contain only digits 

k 

and 1 (for if p contains a digit i > 1 and q contains a digit j > 1, then ix b j = min{i, j} > 1 
will appear somewhere in p x b g). As in the proof of Theorem 12 we may assume that this 
q has the form 11 . . . 1 | 6 for some s with 1 < s < k. Suppose p = ^[Zq pC? with pi 6 {0, 1} 

s 

is a divisor of .11 ... 1. 1 2, so that 

p x 2 (2 fc+1 - r - 1) = 2 k - 1 . (29) 

By Lemma 5, p' := X][=o P»^* ^ s a Dase ^ dismal divisor of ^^ : 

fefc+i-r _ 1 b k _ 1 

" x '^r3^ = i-r- (30) 

Furthermore, (30) still holds if any of the pi that are 1 are changed to any digit in the range 
{1,2,. . . ,b — 1}. Conversely, any base b dismal divisor p' of ^3^- remains a divisor if all the 
nonzero digits in the base b expansion of p' are replaced by l's. So each divisor of 11 ... 1 12 

k 

with t l's corresponds to (b — 1)* divisors of Kzr- Since there are T(k,t) divisors of 2 k — 1 
with t l's, the result follows. □ 

Remarks. 

(i) Table 10 shows the initial values of <4(11 ... 1 1&). 

k 

19 



fc\6 


2 


3 


4 


5 


6 


7 


8 


9 


10 


1 


1 


2 


3 


4 


5 


6 


7 


8 


9 


2 


2 


6 


12 


20 


30 


42 


56 


72 


90 


3 


3 


11 


39 


84 


155 


258 


399 


584 


819 


4 


5 


34 


129 


356 


805 


1590 


2849 


4744 


7461 


5 


8 


82 


426 


1508 


4180 


9798 


20342 


38536 


67968 


6 


14 


206 


1434 


6452 


21830 


60594 


145586 


313544 


619902 


7 


24 


526 


4890 


27828 


114580 


375954 


1044246 


2555080 


5660208 



Table 10: Table of d&(ll . . . 1|&) (with k l's), the number of base b dismal divisors of ^— p 
(rows: A002378, A027444, A186636; columns A079500, A186523). 



(ii) Theorem 13 reduces to Theorem 12 in the case 6 = 2. 
(iii) From (27) and (28) we have 

bk ~ 1 < d b (ll^A\ b ) < (b-l)b k ^ 

k 



k 



(31) 



For b = 2 we also have the asymptotic estimate (23). 



We now study the runners-up in the binary case (among odd numbers of length greater 
than 5), namely the numbers 111 . . . 101 12 and 101 . . . 1 1 1 12- The simplest way to state the 
result is to give the generating function. 



Theorem 14. 




00 
J2d2(2 k -3)z k ~- 

k=3 


Z 3 y, (1 - Z) 2 Z 1 

\-z ^ 1 - 2z + z 1 - 1 - z l + z l + 2 
1=3 




= z + 2z 3 + 2z A + 2z 5 + 4z 6 + 6z 7 + 10z 8 + • • • 



(A188288) (32) 

We will deduce Theorem 14 from Theorem 17 below. 

Suppose p x 2 q = 2 k — 3, k > 3, where len 2 (p) = h, len 2 (g) = I, with h + I = k + 1. In 
order to find all choices for p, we note that the binary expansions of p and q must end with 
. . . 01, and that, as in the proofs of Theorems 12 and 13, we may assume that q = 2 l — 3. 
Our approach is to fix I and allow h to vary. Let M^ denote the number of binary numbers 
p with len 2 (p) = h such that 



px 2 111... 101 1 2 




(33) 



Suppose the binary expansion of p is 



li„Vi • • • x 3 x 2 x 1 01\ 2 , 

where v := h — 3 and the Xj are or 1. The long multiplication tableau for (33) implies that 
the Xi must satisfy certain Boolean equations (remember that + 2 is the logical OR; in what 
follows we will write + rather than + 2 ). For example, the tableau for / = 4 and h = 9, 
v — 6 is shown in Figure 2. 
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1 X6X5X4X3X2X1O 1 

x 2 110 1 

1 X 6 X 5 X4X 3 X2XiO 1 

1 X6X5X4X3X2X10 1 

1 XQX5X4X3X2X10 1 

1111111111 o~T 

Fig. 2. 



By reading down the columns, we obtain the equations 

1. 
1- 
1. 

1. 



X1 + X3 
X\ + £2 + ^4 
X 2 + £3 + £5 
£3 + £4 + Xq 

X 5 + X 6 



1. 



r(4) 



There are Mg = 29 solutions (xi, . . . , a^) to these equations. Table 11 shows the initial 
values of M^ , as found by computer. 



h\l 


1 2 


3 


4 


5 


6 


7 


8 


1 


1 


1 


1 


1 


1 


1 


1 


2 























3 


1 





1 


1 


1 


1 


1 


4 


1 


1 


1 


2 


2 


2 


2 


5 


1 


2 


3 


3 


4 


4 


4 


6 


1 


2 


5 


7 


7 


8 


8 


7 


1 


3 


9 


13 


15 


15 


16 


8 


1 


6 


16 


21 


29 


31 


31 


9 


1 


10 


29 


17 


56 


61 


63 


10 


1 


15 


53 


89 


110 


120 


125 


11 


1 


24 


96 


170 


216 


238 


248 


12 


1 


40 


174 


326 


422 


471 


494 



r(0 



Table 11: Table of M^> (columns 3 and 4 are A070550 and A188223). 
Inspection of the table suggests that the Z-th column satisfies the recurrence 



M, 



(0 = M£. 1 + Mg. a 



1VI h-l+2 + 1VI h-l ^ 1VI h-l-l ) 



(34) 



for / > 3. This will be established in Corollary 18. 

We consider the cases h < l + l and h > 1+2 separately. For h < l+l, it is straightforward 
to show the following: 



M. 



(i) 



0, M, 



(i) 



1 (/i ^ 2) , m; 



(2) 
h 



0. 



(35) 
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and, for I > 3, h < I + 1, 



M, 



V) 



'1, 




if /i = 1, 


o, 




if /i = 2, 


2 h ~ 3 , 




if 3 < /i < I - 1, 


2 h-3 _ 


-1, 


if h = I or / + 1. 



(36) 



This accounts for the entries in Table 11 that are on or above the line h — I = 1. 

We now consider the case 3<l<h — 2 = v + l. The multiplication tableau leads to two 
special equations, 

%i — 1, if / = 3, or 

^x\ + x 2 + ■ ■ ■ +xis + xi-± = 1, if / > 4, 

and 

x v-l+3 "*" x v-l+4 + • • • + X v _\ + X v = 1 , 

together with a family ofv — l + 1 further equations, which, if / = 3, are 



^1 + ^3 = X2+ £4 



X*v— 3 ■ %v— 1 



X v —2 ■" 3-D 



(37) 
(38) 
(39) 



or, if I > 4, are 



£l + £ 2 + X 3 + • • • + X1-2 +Xi = 1 
X2+X3 + X4+ ■ ■ ■ + X1-1 + Xi +1 = 1 



Xv-l+1 + 37jj— Z+2 + ^D-i+3 + • • • + X v -2 + X,; 



1. 



(40) 



The two special equations (37) and (38) involve variables with both low and high indices, 
which makes induction difficult. We therefore define a simpler system of Boolean equations 
in which the special constraints apply only to the high-indexed variables. 

For I > 3 and n > 1, let Dn denote the number of binary vectors £1X2 • • • x n of length n 
that end with x n = 1, do not contain any substring 

.00 . . . 000 or 00 . . . 010 



and do not end with 



Equivalently, D n is the number of solutions to the Boolean equations 

Xi + X 2 + £3 + ■ ■ • + X;_2 + Xi = 1 , 
X 2 + X 3 + Xi+ ■■■ + Xi_i + X t+1 = 1 , 



X n -i + X n -i + i + X n -i + 2 + • + X n - 3 + X n -i = 1 



(41) 
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and 



%n-l+2 + ^n-i+3 + " " ' + ^n-2 + ^n-1 



./. , 



1. 



(42) 



We also set Dq 



(D 



Theorem 15. Fori > 3, n > 1, there is a one-to-one correspondence between binary vectors 
of length n satisfying the D n ' equations and compositions of n into parts taken from the set 



{1,2,. ..,1-2,1,1 + 1}. 



(43) 



Proof. We exhibit a bijection between the two sets. Let c be a composition n = Ci+C2+- • -+c r 
into parts from (43). For 1 < % < I — 2, let ip(i) = 00 ... 01, of length i and ending with a 
single 1, let tp(l) = 00 . . . 011, of length I, let ${1 + 1) = 00 . . . 0011, of length I + 1, and let 



V>(c) = ip(ci) tp(c 2 ) ■ ■ ■ ip{c r ) 



(44) 



a binary vector of length n. Note that the ip(i) for 1 < i < I — 2 contain runs of at most I — 3 
zeros. Runs of / — 2 or / — 1 zeros in ip(c) are therefore followed by two ones. So conditions 
(41) and (42) are satisfied. Conversely, given a binary vector satisfying the Dn equations, 
we can decompose it into substrings ip{i) by reading it from left to right. □ 



n\l 


3 


4 


5 


6 





1 


1 


1 


1 


1 


1 


1 


1 


1 


2 


1 


2 


2 


2 


3 


2 


3 


4 


4 


4 


4 


6 


7 


8 


5 


6 


11 


14 


15 


6 


9 


20 


27 


30 


7 


15 


36 


51 


59 


8 


25 


65 


98 


115 



Table 12: Table of D n l) (the columns are A006498, A079976, A079968, A189101). 

The generating function for the D n follows immediately from the theorem: 
Corollary 16. For I > 3, the numbers D n have generating function 



X>(0 ( 



£^ n = r 



n=0 



'z + Z 2 H h ^~ 2 + z' + z' +1 ) 

1-Z 



l-2z + z 1 - 1 -z l + z l + 2 



(45) 
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Table 12 shows the initial values of D n , computed using the generating function. The 
I = 4 and I = 5 columns are in [17] as entries A079976 and A079968, taken from a paper by D. 
H. Lehmer on enumerating permutations (tti, . . . , ir n ) with restrictions on the displacements 
7Tj — i ([14]; see also [11]). 

We now express the numbers M^ in terms of the D n . We consider the I values xi, . . . , x\ 
in a solution to the M h equations, and classify them according to the number of leading 
zeros. There are just / — 1 possibilities, as shown in Table 13, and in each case the M^ equa- 
tions reduce to an instance of the D n equations. For example, if X\ — 1 the M^ equations 

1, equations (37), 



reduce to an instance of the D h _ 3 equations. (E.g., if I = 5 and we set X\ 



(38), (40) become x v _ 2 + %v-i + x v = I, x 2 + x 3 + x^ + x 6 = 1, . . ., x„_ 4 + x„_ 3 + x„_ 2 + x v 
which, if we subtract 1 from each subscript, are the equations for D v , that is, D h _ 3 .) 



Setting, in Mj /} , 


leads to 


Xi — 1 




D (t) 


X\ = 0,X2 = 1 




D (i) 


xi = x 2 = 0,x 3 = 1 




D (i) 


X\ = ■ ■ ■ = X[_4 = 0, X/_3 


= 1 




Xi = ■■■ = X^ 3 = 0,X Z _ 2 


= Xj-1 = 1 


D (l) 


Xi = ••• = Xi_ 2 = 0,Zj_i 


= Xi = 1 


D (l) 



Table 13: Expressing M^ in terms of D n . 



We have therefore shown that for I > 3, h > 3, 



M; 



(0 



,(0 



,(0 






,(0 



,(0 



D v - 1 4- ■ ■ ■ 4- D { 
u h-h ^ ^ u h-l+\ 



D 



(0 

h-l-l 



D 



(0 

h-l-2 ■ 



(46) 



Table 14 illustrates (46) in the case I = 4. 



h 


Mf 


D {i) 


p(4) 
^-5 


^(4) 


3 


1 


1 


— 


— 


4 


1 


1 


— 


— 


5 


3 


2 


1 


— 


6 


5 


3 


1 


1 


7 


9 


6 


2 


1 


8 


16 


11 


3 


2 


9 


29 


20 


6 


3 



Table 14: Illustrating M, 



(i) 



)( 4 ) 



(4) 



D£U + d%u + D 



)(4) 



From (46) and Corollary 16, and taking into account the values of M^ for h < 3, we 
obtain: 
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Theorem 17. For I > 3, 



M {l \z) := J2 M h )zh = z + (z 3 + z 4 + ■ ■ ■ + z 1 - 1 + z l+1 + z l+2 )V (l \z) 
h=i 

z(l-z) 



1 - + z 2 H h z'~ 2 + z* + ^ z+1 ) 

(47) 



z(i-z) 2 



1 - 2z + z 1 - 1 -z l + z l + 2 

It is now straightforward to obtain the recurrence for M^ from the generating function 
in the second line of the display. We omit the proof. 

Corollary 18. For I > 3, M h satisfies the recurrence (34) with initial conditions (35), 
(36). 

We can now give the proof of Theorem 14. From the definition of M^ , we have 

k 

d 2 (2 k -3) = Y, M i 



k 

l k-l+l ■ 

1=1 



That is, d2(2 k — 3) is the sum of the coefficient of z k in Ai^(z), the coefficient of z k ~ l in 
MS 2 \z)i . . ., and the coefficient of z 1 in MS k '{z). In other words, c?2(2 fc — 3) is the coefficient 
of z k in 

M {1 \z) + zM {2 \z) + z 2 M {? \z) + ■■■ + z k ~ l M {k \z) , 

and now (32) follows from A4^(z) = z + z 3 /(l — z), A^^(^) = 0, and Theorem 17. This 
completes the proof of Theorem 14. 

Remark. The I = 3 column of the M h table (Table 11) is an interesting sequence in its own 
right. 4 To analyze it directly, first consider the system of simultaneous Boolean equations 

xi+x 2 = X2 + X3+ ■ ■ ■ + x n -i +x n = 1 , (48) 

for n > 2, involving a chain of linked pairs of variables. An easy induction shows that the 
number of solutions is the Fibonacci number F n+ 2 (cf. A000045). 5 Second, the equations for 
M^ , (38) and (39), break up into two disjoint chains like (48), and we find that 

M 0) = 1^-2)72^/2, if n is even, 

1 F(n-3)/2F(n+l)/2, ^ U is odd. 



4 It is entry A070550 in [17], which contains a comment by Ed Pcgg, Jr., that it arises in the analysis of 
Pcnncy's game. 

5 This result could also be obtained by the Goulden- Jackson cluster method, as implemented by Noonan 
and Zeilberger [8], [16]. 
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From (49) we can derive the recurrence M^ = M^_ x + M^_ 3 + M^_ 4 and the generating 
function 

*(!-*) 

1 — z — z 3 — z 4 

in agreement with (34) and (47). 



M^\z) = , y J> ^ , (50) 



The final result in this section will describe the asymptotic behavior of the sequence 
G?2(2 fc — 3). When investigating Conjectures 3 and 4, we observed that among all numbers 
n with k binary digits, the number 2 k — 1 was the clear winner, with values close to the 
estimate (23). The runners- up, a long way behind, were 2 fc — 3 and 2 k — 2 k ~ 2 — 1 (and 
sometimes other values of n), all with the same number of dismal divisors, for which the 
number of dismal divisors appeared to be converging to one-fifth of the number of divisors 
of the winner, or in other words it appeared that 

4(2 fc -3) 1 

4(2^1) "> 5> as ^°°- ( 51 ) 

We will now establish this from the generating function (32). 

Theorem 19. 

d 2 (2 k -3) ~ -4— r (1 + 9*), ask-^oo, (52) 

ok log 2 

where Ok is a bounded oscillating function with |G&| < 10 



-5 



Proof. Our proof is modeled on Knopfmacher and Robbins's proof [12] of (23), which uses 
the method of Mellin transforms as presented by Flajolet, Gourdon, and Dumas [5]. We will 
indicate how the Knopfmacher- Robbins proof can be reworded so as to establish (23) and 
(52) simultaneously. 

Knopfmacher and Robbins work, not with (22), but with 

/(*) == V {1 - Z)Z \ , (53) 

J w 4^ 1 - 2z + z l y ' 

which is the generating function for the number of compositions of n into parts of which the 
first is strictly greater than all the other parts (A007059). Equations (22) and (53) basically 
differ just by a factor of z. Then [12] shows that the coefficient of z 11 in f(z) is 

')n— 1 

1 + 6) , (54) 



nlog2 



for some small oscillating function O, which implies (23). 

So as to have a function with the same form as (53), we consider, not (32), but 

*— ' 1 — 2z + z l — z l+l + z l+6 
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and we will show that the coefficient of z n is 

2 n+l 



5n log 2 



1 + 0) , (56) 



for some (different) small oscillating function 6, which implies (52). We can change the 
lower index of summation in (55) from 2 to 1, since the I = 1 term is the generating function 
for the Padovan sequence (A000931), which grows at a much slower rate than (56). 

In what follows, we simply record how the expressions in Knopfmacher and Robbins's 
proof [12] need to be modified so as to apply simultaneously to (53), which we refer to as case 
I, and (55), which we call case II. We follow Knopfmacher and Robbins's notation, except 
that we use j and m as local variables, rather than k, to avoid confusion with the k in the 
statement of the theorem. Some typographical errors in [12] have been silently corrected. 

Let pj denote the smallest root of the denominator of the j-th summand in f(z) that lies 
between and 1. Then 

Pi = l(l + TZ-* + 0(j2-*)), (57) 

where r = 1 (case I) or 5/8 (case II). 

Let q n j denote the coefficient of z n in the j-th summand in f(z). Then 

q n ,j « T~^ (1 - J)" 

ps T~ 3 ~ e e~ Tn/23 , (58) 

where e = 1 in case I or 2 in case II. Next, f n , the coefficient of z n in f(z), is 

2 "~ e ( ]T2^e- rn/2J +o(l) j . 



Let 



oo 



g{x) := J2 2 ~ 3e ~ TX/23 



i=2 
The Mellin transform of g(x) is 

1 o 2 ^ 1 ) 
9*{*) ■= ^7-T^T r ( s )' 0<»( a )<l. 

To compute f n we use (following Knopfmacher and Robbins) the Mellin inversion formula 

i ^l/2+ioo 

9(x) = 7T~ x~ s g*(s)ds. 

2m Jl/2-ioo 

Now g*(s)x~ s has a simple pole at s = 1 + Xm, for each m G Z, where Xm = 2irim/ log 2, 

with residue 

1 r(l + x m )e- 27rimlog2X 

x log 2 r l+Xm 
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After combining the contributions from all the poles, we have 

1 T(l + 2nim/\og2)e~ 2mmlo ^ n 



In 



£ 



m=— oo 



H log 2 T l+2mm/ log 2 

The term for m = dominates, and we obtain the desired results (54) and (56). 



(59) 
□ 



7 The sum of dismal divisors 

We briefly discuss the dismal sum-of-divisors function <7b(n) (see A188548, A190632, and 
A087416 for bases 2, 3, and 10). 

Theorem 20. In any base b>2, z/Ten^n) = k, then 



and <Jb{n) — n if and only if n 



n < cr b (n) < b k 
b — 1 (mod b). 



(60) 



Proof. The first assertion follows because n divides itself, and no divisor has length greater 
than k. If n ^ b — 1 (mod b) then since b — 1 is the multiplicative unit, <7&(ra) ^ n. Suppose 
n = b — 1 (mod b) and p is a divisor of n, with say pX b q = n. Both p and g must end with 
(3. From the long multiplication tableau, p -C;, n, so p+b n = w, and therefore cr&(n) = n. □ 

In ordinary arithmetic, a number n is perfect if its sum of divisors is 2n. In dismal 
arithmetic, n+^n = n. So the second part of the theorem might, by a stretch, be interpreted 
as saying that the numbers congruent to b— 1 (mod b) are the base b perfect dismal numbers. 

In base 2, then, 02(71) = n if and only if n is odd. Examination of the data shows that, 
if n is even, with len 2 (n) = k, often 02(71) takes its maximal value, 2 fc — 1. Table 15 shows 
the first few exceptions, which are characterized in the next theorem. 



n 


^(n) 


n 


02 («) 


10010 


11011 


1001010 


1101111 


100010 


110011 


1001110 


1101111 


100110 


110111 


1010010 


1111011 


110010 


111011 


1100010 


1110011 


1000010 


1100011 


1100110 


1110111 


1000100 


1110111 


1110010 


1111011 


1000110 


1100111 







Table 15: Even numbers n such that o~2{n) is not of the form 11 . . . 1| 2 (A190149-A190151). 
Both n and 02(71) are written in base 2. 

Theorem 21. Suppose n = 2 r m with r > 1, m odd, and len 2 (n) = k. Then 

2 (n) = 2 fc -l = 11_^_1| 2 

k 

unless the binary expansion of m contains a run of more than r consecutive zeros. 
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Proof. Since m is odd, 02 (m) = m. Therefore 

cr 2 (n) = m| 2 +2 ttz.0| 2 +2 ^00| 2 +2 •■■ +2 m QO . . . 0J 2 , 

r 

and any string 00 ... 1 in m will become 11 ... 1 1 in 02(71) unless i exceeds r. D 

i i 

The first entry in Table 15 is explained by the fact that n = 100 10 1 2, r = 1, m = 1001 1 2, 
and m contains a run of two zeros. 

We also considered two other possible definitions of perfect numbers: (i) n is perfect in 
base b > 3 if the dismal sum of the dismal divisors of n is equal to 2 x b n. We leave it to 
the reader to verify that for this to happen, b must be 3, and then n is perfect if and only 
if n = 2 mod 4. (ii) n is perfect in base b > 2 if the dismal sum of the dismal divisors of n 
different from n is equal to n. But here n cannot be b — 1, so b — 1 is a divisor, and n ends 
with 6—1. This implies that n has no divisors of length len^n) except n itself, so the sum 
cannot equal n, and therefore no such n exists. 

We end this section with a conjecture (see A186442): 

Conjecture 5. For all n > 1, dxo(n) < cr 10 (n). 

8 Dismal partitions 

Since n +t,n = n, it only makes sense to consider partitions into distinct parts (otherwise 
every number has infinitely many different partitions). We define p&(n) to be the number of 
ways of writing 

n = m 1 + b m 2 + b ■■■ + b mi, (61) 

for some I > 1 and distinct positive integers rrii, without regard to the order of summation. 
We set pb(0) = 1 by convention. 

For example, Ps(7) = p 3 (21\ 3 ) = 22, since (working in base 3) 21 is equal to 21 +3 any 
subset of {20, 11, 10, 1} (16 solutions), 20 +3 11 +3 any subset of {10, 1} (4 solutions), and 
20 +3 1 +3 any subset of {10} (2 solutions), for a total of 22 solutions. 
Remarks. 

(i) Permuting the digits of n does not change pb{n). 

(ii) Any zero digits in n can be ignored. If n' is the base b number obtained by dropping 
any ?Vs that are zero, Pb(n') = Pb(n). 

(iii) Although we will not make any use of it, there is a generating function for the pb{n) 
analogous to that for the classical case. If we interpret z m x b z n to mean z m bn , then we 
have the formal power series 

l+Pb(l)z + p b (2)z 2 +p b (3)z 3 + --- = (l + z)x b (l + z 2 )* b (l + z 3 )x b ■■■ . 

(iv) The sequences P2{n) and pio(n) form entries A054244 and A087079 in [17], con- 
tributed by the second author in 2000 and 2003, respectively. 

In the remainder of this section we index the digits of n by {1,2, . . . ,n}, in order to 
simplify the discussion of subsets of these indices. 
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Theorem 22. If n — n\n 2 ■ ■ -n k \ 2 , ^i G {0, 1}, and the binary weight of n is w, then P2(n) 
is equal to the number of set-covers of a labeled w-set by nonempty sets (cf. A003465, [4, 
p. 165]), that is, 



»(»> = \ B-ir-Qa 2 ' 



(62) 



Proof. From Remark (ii), p2{n) = P2QI • • • 1 h)- There is an obvious one-to-one corre- 

w 

spondence between collections of distinct nonempty subsets of {l,...,w} whose union is 
{1, . . . , w} and sets of distinct nonzero binary vectors whose dismal sum is 11 ... 1 [2- □ 

w 

Theorem 23. If n = n±n 2 . . . n k \b, < n* < 6 — 1, then 

p b (n) = \ Y. (-1) |5| 2^K+^, (63) 

sc{i,...,k} 

where e^ = if i G S , q = 1 if i ^ S . 

Proof. The set Q n of x <C{, n is a partially ordered set (with respect to the operator <C?,) 
with Mobius function given by [21, §3.8.4] 

f(-l)Ei(2/i-^i) if y . - x . = or 1 for alii, 

fx(x 1 ...x k \ b , y 1 ...y k \ b ) = < (64) 

10, otherwise. 

Every subset of Q n \ {0} has dismal sum equal to some number <C;, n, so we have 

J2 Pb{x) = 2^^+!)-! . 

From the Mobius inversion formula [21, §3.7.1] we get 



I y^ (_i)i^i2n ies ^n^sK'+ i ) 



Pb{n) 

sc{i,...,k} 

which implies (63). D 

Theorem 23 reduces to Theorem 22 if all n« are or 1. 

Corollary 24. Suppose n = n\n 2 . . .nk\b, and let y be the ordinary product n\ x n 2 x • • • x n k . 
Then pb{n) is divisible (in ordinary arithmetic!) by 2 y ~ 1 . 

Corollary 25. For a single-digit number n = n\\ b , Pb{n) = 2™ 1 " 1 . For a two-digit number 

n = n l n 2 \b, 

p b (n) = 2 (ni+1){n2+1) ~ 2 - 2 nin2 ~ 1 (2 ni -l)(2 n2 -l). (65) 

Eq. (65) was found by Wasserman [17, entry A087079]. 

It follows from the above discussion that in any base b, the only numbers n such that 
Pb{n) = 1 are 0\ b , l\ b , 10 1&, 100 |j,, 1000|&, . . ., and that all other numbers n have the property 
that the dismal sum of the numbers x <C& n is n. These two classes might be called "additive 
primes" and "additive perfect numbers." 
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9 Conclusion and future explorations 

We have attempted to show that dismal arithmetic, despite its simple definition, is worth 
studying for the interesting problems that arise. We have left many questions unanswered: 
the "prime number theorem" of Conjecture 1, the questions about the numbers of divisors 
stated in Conjectures 2-4 (in particular, is it true that 11 . . . 1 |io has more divisors than 
any other base 10 number with the same number of digits?), the base 10 dismal analog of 
d{n) < a(n) (Conjecture 5), and the two questions about dismal squares at the end of §4 — in 
particular, is there a recurrence for the sequence (19)? There are numerous other questions 
that we have not investigated (for example, if x <Cf, y, what can be said about x Xi,y7). 

We have made no mention of the complexity of deciding if a number is a dismal prime, 
or of finding dismal factorizations. In base 2 such questions reduce to solving a set of 
simultaneous quadratic Boolean equations, where a typical equation might be 

(^o*2 2/4) +2(3:1X2^/3) +2 ••• +2(3:4*2 yo) = 0(orl). 

This becomes a question about the satisfiability of a complicated Boolean expression, and is 
likely to be hard to solve in general [7]. 

While we have focused on the cases b = 2 and b = 10, it would be nice to better 
understand the qualitative differences across a wider range of bases. For example, while 
b = 2 is a kind of Boolean arithmetic, does 6 = 3 correspond to a three-valued logic? Do 
odd b behave differently from even W. More generally, what other interesting mathematical 
structures might be modeled by dismal arithmetic? 
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